Stabilization/destabilization of cell membranes by multivalent ions: Implications for 

membrane fusion and division 
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We propose a mechanism for the stabilization/destabilization of cell membranes by multivalent 
ions with an emphasis on its implications for the division and fusion of cells. We find that multivalent 
cations preferentially adsorbed onto a membrane dramatically changes the membrane stability. They 
not only reduce the surface charge density of the membrane but also induce a repulsive barrier to 
pore growth. While both of these effects lead to enhanced membrane stability against vesiculation 
and pore growth, the repulsive barrier arises from correlated fluctuations of the adsorbed cations 
and favors closure of a pore. Finally, the addition of a small amount of multivalent anions can 
reverse the membrane stabilization, providing an effective way to regulate membrane stability. 
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Lipid bilayer membranes are resistant to rupture, pri- 
marily serving as a barrier to the leakage of the cell's 
contents, while also being dynamic structures that un- 
dergo various topological transitions. The capability 
of living cells to regulate the stability of their bound- 
ing membranes is crucial to their maintenance and re- 
production [Q. Membrane stability against rupture 
changes most dramatically during cell division and fu- 
sion. The precise mechanism for achieving this complex 
task in living cells is complicated by various membrane- 
associated/bounded proteins M and is not yet clear. Nu- 
merous studies, however, suggest that membrane sta- 
bility is influenced by several factors such as the ionic 
strength, external fields, and thermal fluctuations [p|-^2|. 
For example, red blood cells can be converted into vesi- 
cles by osmotic lysis in a solution of low ionic strength 
lacking multivalent cations The presence of di- 

valent cations, however, prevents this vesiculation [^-^). 
In fact, a number of experiments j4[j|] have unambigu- 
ously demonstrated that the stability of red cell mem- 
branes against vesiculation can be greatly enhanced by 
multivalent cations. Despite this, a consistent theoretical 
description of this phenomenon has so far been lacking. 

The strong valency dependency of membrane stabil- 
ity motivated this work. Not only can osmotic lysis 
lead to vesiculation, but it can also create large pores in 
the cell membranes that subsequently contract to a size 
that is controlled by the ionic strength. Pore closure can 
be stimulated by cations, and remarkably the rate of pore 
closure strongly depends on the valency of cations ||; 
Ca 2+ is roughly 60 times as potent on a molar basis as 
Na + . The potency of divalent cations, which essentially 
prevents vesiculation, was first demonstrated experimen- 
tally three decades ago Q , but its has yet to be examined 
theoretically. Here we propose a theoretical mechanism 
to explain this phenomenon. We find that multivalent 
counterions adsorbed onto charged membranes dramat- 
ically enhance the membrane stability through two ef- 
fects. First, they can significantly reduce the strength of 
the electrostatic repulsion between backbone charges on 



the membrane, which enhances the membrane stability. 
Second, they induce a repulsive barrier to pore growth. 
The repulsive barrier originates from the correlated fluc- 
tuations of adsorbed counterions and favors closure of a 
pore, further stabilizing the membranes against vesicu- 
lation and pore growth. Upon adding a small content 
of multivalent anions, the adsorbed cations are released 
into solution, thus reversing the membrane stabilization. 
Adsorption/desorption of multivalent cations provide an 
effective way to regulate the membrane stability. 

The model we consider here is a thin flat mem- 



brane 1 13 1 in the xy-plane, in the presence of a mono- 
valent (1:1) salt such as NaCl, and in the presence or 
absence of Z + -valent (Z + : 1) salts, such as CaCl2. Each 
side of the membrane is assumed to be negatively charged 
with constant charge density — eoo and attracts ions of 
the opposite charge, as schematically shown in Fig. [l](a). 
For simplicity, we consider the case of a single circular 
pore of radius R, already formed in the membrane by, 
for example, osmotic stress. The stability of the mem- 
brane against rupture can be quantified in terms of a 
line tension, i.e., the energetic penalty for creating a 
pore per unit length. The electrostatic repulsion between 
charges on the membrane favors pore formation J9|, but 
the hydrophobic effect tends to close the pore. If mobile 
ions are treated as screening objects that simply reduce 
the electrostatic repulsion between charges on the mem- 
brane via Debye screening, then the electrostatic contri- 
bution to the line tension ^dh can be estimated using 

Debye-Huckel (DH) theory ft: lDH e^a^K^e^R, 

ifR< k- 1 and lDH _^JX.-2^-i it 
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if R > k~ , where 

k^ 1 is the screening length and e is the dielectric constant 
of the solvent. Obviously, the electrostatic repulsion fa- 
vors creation and expansion of a pore, i.e., jdh < 0. 

Charged membranes are, however, capable of ad- 
sorbing counterions of the opposite charge (See 
Fig. 0(a)) p^ , [l5|l, reducing the surface charge density of 
the membrane |14|]. The magnitude of the reduced, renor- 
malized charge density can be estimated by equating the 
chemical potentials of the "free" and "condensed" coun- 
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terions, i.e., those adsorbed onto the membrane surface. 
In the following descriptions, the subscripts i = 1 and 2 
refer to the monovalent and multivalent counterions, re- 
spectively. If o"i is the number density of condensed coun- 
terions, then the effective (renormalized) surface charge 
density on the membrane is — ea* = — e(<jQ — o\ — Z+o-i)- 
The chemical potential of free counterions is mainly 
associated with the configurational entropy of mixing: 
^/ree ^ ksT hi (maf\ where rii and ai are the concen- 
tration and size of counterions, respectively. On the other 
hand, the chemical potential of the condensed counteri- 
ons arises from electrostatic interactions and the entropic 

penalty for condensation; pf° nd k B T(Z + £ B a*VS) + 

In {(Jiaf) if kT x > y/~S, where S is the area of the mem- 
brane, and [f° nd k B T(Z + l B cy* n- 1 ) + \n{a l af) oth- 
erwise. The equilibrium values of cr, can then be obtained 
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FIG. 1. (a) Schematic view of a charged membrane with 
a pore of a radius 7?. Charged and neutral lipids are de- 
noted by — and 0, respectively, while the adsorbed multiva- 
lent cations are denoted by encircled +'s. A pair of attracting 
lipid and lipid/cation is also shown that stabilizes the mem- 
brane against pore growth, (b) Illustration of multivalent 
anions (encircled — 's) forming pairs with the adsorbed mul- 
tivalent cations, then leaving into solution. Depletion of the 
adsorbed layer destabilizes the membrane. 



The condensed counterions do not simply renormal- 
ize the membrane charge density but also give rise to 
charge fluctuations in the plane of the membrane surface 
that tend to be correlated with each other ]l6| , p^ |. A 
typical attracting pair of a lipid and lipid/cation is illus- 
trated in Fig. [I]. Creation of a pore makes charges at 
the edge less efficiently correlated and is discouraged by 
the charge correlation effects. Computation of the charge 
correlation contribution to the pore free energy is highly 
involved, partly because the charge fluctuation interac- 
tions are not pairwise additive flqj[7|| . This is compli- 
cated by yet another factor: the specific geometry of the 
membrane with a pore. The electrostatic effects at the 
meanfield level, suppressing both adsorption and charge 
fluctuations, has only recently been addressed To 
study the effects of charge fluctuations on the membrane 
stability, we take the continuum limit and incorporate 



the in-plane charge fluctuations at the Gaussian level as 
in previous cases [|l6]jlj]. The resulting charge fluctua- 
tion contribution to the pore free energy, i.e., the change 
in the charge fluctuation free energy by creating a pore, 
is formally given by 
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Here, the matrix Q is defined by the matrix elements 

e -«|x x -x^| 



1 + ^BTccCx 



X' 



(2) 



where xj_ = {x,y), i B — e 2 /ek B T is the Bjerrum length, 
a cc = o\ + Z + <T2, k 2 = &Ttl B I, I is the ionic strength of 
the solution, and Q = lim/j^o Q. Finally Cx^x^ = 1; if 
xj_ and x^ are on the membrane and is otherwise. In 
the case of S > as is the case for red cell experi- 

ments [|]||, J-pore in Eq- (Q) can be calculated without 
making further approximations. This follows from the 
fact that 8 = (Q — Q)Q~ 1 ~ S^ 1 and !F pore can then be 
expanded in powers of 6. In the case S > K , we can 
take the limit S — > oo without introducing any apprecia- 
ble error. In this limit, only the leading term survives in 
the expansion. 

To calculate the free energy in Eq. (|l|) , it proves useful 
to Fourier transform it from xj_ to kj^ . We find 
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o 1 + X cc \/k 2 ± + n 2 \Jr 2 + r' 2 - 2rr' cos 9 
xJoik^r 2 +r> 2 -2rr' cos 9) (3) 



where A cc = l/2ir£ B a cc and Jo (a;) is the zeroth-order 
Bessel function of the first kind. In the case of 
kR 3> 1, the charge-fluctuation line tension, i.e., A7 = 
AJ r pore /2irR, shows two distinct scaling behaviors: A7 ~ 
\~} ln(l/A cc Ac) for small A cc n and A7 ~ \~ 2 k _1 for large 

X CC K. 

The charge fluctuation (A7 > 0) and hydrophobic con- 
tributions (70 > 0) favor closure of the pore, while the 
DH (jdh < 0) or renormalized DH term (7^ < 0) 
tends to expand the pore. To study the membrane stabil- 
ity, we have solved for a* and jtotai = 7o +7* {° t 1dh) + 
A7 simultaneously. We have chosen a\ = a 2 = 2k 
T = 300, and ao = 0.2nm~ 2 , which is in the physio- 
logical range. To suppress the i?-dependency, we have 
assumed that R > Fig. || shows jtotai in units of 

70 = 10~ n J/m ||], as a function of the monovalent coun- 
terion concentration n\. In the absence of multivalent 
counterions (Z = 1), "/total is negative when n\ is in the 
range < n\ < ImM. This implies that the membrane is 
unstable to pore growth as long as n\ is in this range. The 
presence of O.lmM of multivalent counterions (Z = 2 and 
3), however, dramatically enhances the membrane stabil- 
ity. In this case, jtotai is positive for the whole range of 
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n\ and is in the range O.670 < jtotai < O.770. In order to 
enhance the membrane stability up to this level by mono- 
valent counterions, about 5mM concentration would be 
needed. This is approximately 50 times higher than that 
of the divalent counterion concentration. This estimate 
is remarkably consistent with the experimental finding 
that Ca 2+ is roughly 60 times more effective on a molar 
basis than Na + in stimulating pore closure || . 
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FIG. 2. Total line tension, in units of the hydrophobic 
contribution 70 = 10 _11 J/m, as a function of the monova- 
lent salt concentration ni. We have chosen T = 300 and 
cto = 0.2nm -2 . In the absence of multivalent cations (Z — 1), 
there exists a finite range of the monovalent salt concentration 
where the membrane is unstable ("/total < 0) to pore forma- 
tion. The presence of as small a concentration as O.lmM 
of multivalent cations (Z = 2, 3) stabilizes the membrane 
against pore growth for the whole range of monovalent salt 
concentration. The distinction between the monovalent and 
multivalent cases is, however, minor in the DH approach, and 
0.1 mM of multivalent cations only slightly enhances the mem- 
brane stability. 

Our results are striking; the presence of multivalent 
counterions is more crucial to the membrane stabil- 
ity than that of the monovalent salt, though the ionic 
strength is mainly determined by the latter. This can, 
however, be understood in the context of "counterion 
condensation" . The long-ranged electrostatic interac- 
tions allow the membrane to adsorb multivalent counteri- 
ons preferentially even when n 2 < ni. These condensed 
multivalent counterions not only reduce the repulsion 
between backbone charges on the membrane, but they 
also enhance the strength of charge fluctuations — both of 
these effects are more efficient with multivalent cations 
than with monovalent ones. When combined, these two 
effects lead to significantly enhanced membrane stabil- 
ity against pore growth and vesiculation. Also note that 
trivalent counterions are even more efficient in enhanc- 



ing the membrane stability than the divalent counterions. 
As shown in the figure, the dramatic distinction between 
the monovalent and multivalent cases is missing in the 
DH theory. The enhanced membrane stability by multi- 
valent counterions seen in the experiments pLp| can be 
explained only when both the preferential adsorption of 
multivalent counterions and the effects of charge correla- 
tions are properly taken into account. 

Whether a pore grows or closes also depends on the 
height of the barrier as a function of the pore size R. In 
Fig. ||, we have plotted the pore free energy as a func- 
tion of R, in units of fcgT. We have chosen T — 30QK 
and (To = 0.2nm~ 2 . The barrier height is finite in the 
presence of monovalent ions only (Z = 1). In contrast, 
the pore free energy in the presence of O.lmM of multi- 
valent counterions (Z — 2,3) grows indefinitely with R. 
This implies that formation of a large pore is energet- 
ically greatly disfavored, in the presence of multivalent 
counterions. The results are indeed consistent with the 
experimental observation that the presence of O.lmM of 
MgS0 4 stabilized the red blood cell ghosts against vesic- 
ulation Q; a pore originally created by osmotic stress 
will grow into a large one or shut down, depending on 
the buffer quality and the strength of the restoring force 
provided by the spectrin network. Note that the osmotic 
stress will eventually be removed. When Z = 1, the pore 
can grow into a large one, once the barrier is overcome 
by osmotic stress. This will lead to vesiculation if the 
restoring force is outweighed by the repulsion between 
the charged groups on the red cell membranes. In con- 
trast, there is a subsequent barrier to pore growth in the 
case Z = 2,3. This prevents the ghosts from breaking 
into vesicles. In contrast, the DH approach mistakenly 
implies that the barrier height is roughly insensitive to 
the valency of counterions. Thus our results in Fig. |3| 
further support the importance of the counterion valency 
and charge correlations to the membrane stability, con- 
sistent with experiments 

The fact that multivalent cations can be preferentially 
adsorbed onto a charged surface implies that the layer of 
the adsorbed cations can be depleted by multivalent an- 
ions. Imagine an anion of valency Z_ making a pair with 
a cation in the condensed layer and leaving into the solu- 
tion, as illustrated in Fig. |l|(b). Whether this is feasible 
can be tested by calculating the change in the chemical 
potential: A/i ~ -Z + Z-£ B / {a 2+ + a 2 -) + Z + £ B o-* + 
In (712-0-2- / n 2+»2+) ! where the subscript 2 refers to mul- 
tivalent ions and the subscripts + and — refer to the 
cations and anions, respectively. When the valency of 
anions is sufficiently high, this change can be negative. 
This implies that a certain fraction of the cations in the 
layer make pairs with multivalent anions and will return 
to the solution to maintain "chemical equilibrium" . Note 
that a minimum concentration of multivalent anions is 
required for this process to occur. This can be readily 
seen by taking the limit n 2 - — ► and noting that, in this 
limit, the entropic penalty for pairing is too large. Since 
only a very small concentration of multivalent cations 
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is needed to enhance the membrane stability, the pres- 
ence of an equally small content of multivalent anions 
suffices to deplete the condensed layer, effectively reduc- 
ing the membrane stability. Cells contain multivalent 
cations (e.g. Ca 2+ and Mg 2+ ) as well as multivalent an- 
ions (e.g., PO^ - and anionic proteins) and thus could 
regulate their stability in this way. This may provide 
a new insight into the biological phenomena of "break- 
down" and "reassembly" of a nuclear envelope during cell 
division. During cell division, a nuclear envelope disin- 
tegrates into vesicles, which eventually reassemble into 
daughter cells jjj . This is suggestive of the cyclic change 
in the membrane stability: stable — ► unstable — > stable. 
The stabilization/destabilization of membranes by mul- 
tivalent ions could be relevant to this flow of membrane 
stability in a dividing cell. This is, however, complicated 
by various membrane-bounded proteins GJ] and further 
consideration is certainly warranted. 



necessarily lead to membrane fusion and that membrane 
fusion (in the absence of fusion peptides) should be fol- 
lowed by lateral phase separation of the lipid molecules 
into two distinct phases |15| : anionic lipid-poor phases 
and anionic lipid-rich phases "coated" with multivalent 
cations. It is the uncoated phases that undergo a topo- 
logical change and eventually fuse into each other, con- 
sistent with our picture of multivalent cations as efficient 
agents for stabilizing membranes against rupture. 

We have benefited from useful discussions with D. 
Boal, M. Wortis, S. Davies, A. Rutenberg, M. Howard, 
K. Delaney, and G. Tibbits. This work was supported by 
the Natural Sciences and Engineering Research Council 
of Canada. 
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FIG. 3. Pore free energy as a function of the pore size R. 
We have chosen the same parameters as were used to gener- 
ate Fig. ^. The pore free energy is estimated in units of fcsT. 
For the monovalent salt case (Z = 1), the pore free energy 
has a finite barrier. In the presence of O.lmM of multivalent 
counterions (Z = 2, 3), the pore free energy grows indefinitely 
with R. In the DH approach, the barrier height is roughly 
insensitive to the counterion valency Z. 

To summarize, we have presented a working mecha- 
nism for the stabilization/destabilization of membranes 
by multivalent ions. Multivalent cations preferentially 
adsorbed onto a charged membrane not only reduce the 
surface charge density of the membrane but also induce a 
repulsive barrier to pore growth that favors the closure of 
a pore. The main advantage of the membrane stabiliza- 
tion by multivalent cations lies in that this can be easily 
reversed; the addition of a small concentration of multiva- 
lent anions can reverse the membrane stabilization. Our 
results are also consistent with the experimental findings 
that membrane adhesion by multivalent cations does not 



[1 
[2 

[3. 

[4 

[5, 

[6 

[<: 

[8 
[9 

[10 

[11 
[12 
[13 



[14 
[15 

[ie; 
[17; 
[is; 

[19 



B. Albert et al., Molecular Biology of the Cell (Garland 
Publishing, Inc., New York, 1989). 

E.A. Neumann, A.E. Sowers, and C.A. Jordan, Electro- 
poration and Electrojusion in cell biology (Plenum Press, 
New York, 1989). 

D.C. Chang, B.M. Chassy, and J. A. Saunders, Guide to 
Electroporation and Electrofusion (Academic Press, New 
York, 1992). 

T. Steck et al., Science 168, 255 (1970). 

M.L. Lieber and T.L. Steck, J. Biol. Chem. 257, 11651 

(1982); J. Biol. Chem. 257, 11660 (1982). 

V.L. Lew, A. Hockaday, C.J. Freeman, and R.M. 

Bookchin, The Journal of Cell Biology 106, 1893 (1988). 

M. Winterhalter and W. Helfrich, Phys. Rev. A 36, 5874 

(1987). 

D. V. Zhelev and D. Needham, Biochimica et Biophysica 
Acta 1147, 89 (1993). 

M.D. Betterton and MP. Brenner, Phys. Rev. Lett. 82, 
1598 (1999). 

JD. Moroz and P. Nelson, Biophys. J. 72, 2211 (1997). 
W. Sung and P.J. Park, Biophys. J. 73, 1797 (1997). 
J.C Shillcock and D.H. Boal, Biophys. J. 71, 317 (1997). 
More realistically, a cell membrane is curved with unbal- 
anced ionic strength between the inside and outside of 
the membrane. We consider the stage where the osmotic 
stress has been removed via ion transport through the 
pore. 

S. Alexander et al. J. Chem. Phys. 80, 5776 (1984). 
J. Israelachvili, Intermolecular and Surface Forces, Chap- 
ter 12 (Acamedic Press, San Diego, 1991). 
B.-Y. Ha and A.J. Liu, Phys. Rev. Lett. 81, 1011, (1998); 
ibid 79, 1289 (1997) Phys. Rev. E 60, 803 (1999). 
R. Podgornik and V. A. Parsegian, Phys. Rev. Lett. 80, 
1560 (1998). 

P.A. Pincus and S.A. Safran, Europhys. Lett. 42, 103 
(1998). 

J.N. Israelachvili and Patricia M. McGuiggan, Science 
241, 765 (1988); D.E. Leckband, C.A. Helm, and J. Is- 
raelachvili, Biochemistry 32, 1127 (1993). 



4 



